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I 1 Abstract 
Oh" 

' We study the exclusive production of heavy flavors at central rapidities in hadron-hadron 

collisions within the kr factorisation formalism. Since this involves regions of small Bjorken x 
in the unintegrated gluon densities, we include the next-to-leading order BFKL contributions 
working directly in transverse momentum representation. Our results are presented in a form 
suitable for Monte Carlo implementation. 

CN 

1 Introduction 

<N 

Scattering processes with at least one hard scale are typically well described using perturbative QCD 
in the framework of collinear factorisation. In this approach cross sections are written as a convolution 
Q\ [ of a purely perturbative partonic cross section with non-perturbative parton distribution functions. 

The latter follow the DGLAP evolution which describes their dependence on the hard perturbative 
scale. When the center of mass energy is very large compared to the perturbative hard scale, or a final 
state is fixed such that there are large rapidity differences among the emitted particles, an alternative 
high energy factorisation based on the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution equation 
applies [I]. Here the hard subprocess is convoluted with the hadron structure using unintegrated 
gluon densities which include their kr dependence in the small Bjorken x limit. This can be seen 
as the small x limit of kx factorisation. In this case enhanced logarithmic in x contributions are 
resummed. 

In the present letter we propose to take the exclusive production of heavy quark-antiquark pairs 
in the early data at the LHC as a test of this formalism and the use of unintegrated gluon densities. 
Large masses such as those of bottom or top quarks allow for a perturbative treatment. In the case of 
top quark pairs their masses are so large that the typical proved values of Bjorken x are not that small. 
In this case it is known that cross sections receive significant corrections from threshold logarithms 
(see [2j[3l|4l[5j[6l[7l[8l[9l [10] for recent results in this direction) . The bottom quarks are lighter and 
therefore test regions of smaller values of x where the corresponding resummations find their natural 
environment. Previous investigations of heavy quark production similar to our present calculation 
where presented in [11]. What we will show in this letter is an alternative approach which operates 
with NLO unintegrated gluon densities in transverse momentum space, does not involve the use of 
anomalous dimensions, treats the kinematics of the quark-antiquark pair exclusively and is readily 
suitable for a Monte Carlo analysis which we will present elsewhere. Other works which we found of 
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interest in the field of inclusive heavy flavor production are Refs. (TTJ [12l [T3l [Ml |T5l [16j ttZ] • For the 
production of bottom pairs in particular we highlight Rcfs. [HI [19], where the reported agreement 
with experimental data at the Tevatron ranges from reasonable [18] to very good [19] , 

The fully exclusive study that we propose could be very useful at the LHC since it allows for 
the precise determination of the x values at which the unintegrated gluon densities are probed. This 
provides a good control on the accuracy of the approximations that we use in our calculation. At the 
LHC the dominant production process for both top and bottom pairs is given by gluon-gluon fusion. 
However, as we already pointed out, only the bottom pair production occurs at relatively small x 
providing the correct kinematics to apply high energy kx factorisation. Top pair production, on the 
other hand, occurs at relatively large values of x due to the large top mass. Studies of its exclu- 
sive production would certainly require the matching of the present calculation with renormalization 
group evolution and can be therefore considered as a test of the capability to extend our high energy 
factorisation towards the region of large x. For our predictions we incorporate the NLO corrections 
to the BFKL evolution kernel [20]. A related study, devoted to the exclusive central production of 
jets in hadron-hadron collisions in kx factorisation, was presented in Ref. |21j . 

After this brief Introduction, in Section 2, we present the general structure of the kx factorised 
differential cross-section and calculate its different elements. In Section 3 we discuss the unintegrated 
gluon density in kx space and its iterative structure. Finally, we write our Conclusions in Section 4. 



2 The kx factorised differential cross-section at NLO 

To describe the differential cross-section for the exclusive production of a pair of heavy quarks within 
kx factorisation it is convenient to introduce a Sudakov basis. To this end we define the light-like 
momenta p\ and P2 which coincide in the s — > oo limit with the momenta of the incoming protons pa 
and p b ■ 

mp mp 

Pi = Pa Pb, P2 = Pb Pa, (1) 

s s 

with s = (pa + Pb) 2 being the squared center of mass energy of the hadronic process. With these 
definitions, we can then work with the usual Sudakov decomposition of a general four momentum, i.e. 

k = a Pl +f3p 2 + k ± . (2) 




Figure 1: Central production of two heavy quarks in kx factorisation 



The notation for the relevant momenta in the partonic hard subprocess is given in Fig. [TJ In the 
BFKL formalism i-channel gluons carry a modified propagator which reggeises them. This propagator 
is associated to the momenta q\ and qi in Fig. [TJ These simplify in the high energy limit and can be 
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written as 

qi = xipi+q^x, 92 = x 2 p2 + Q2,±- (3) 

On the other hand, the momenta of the produced heavy quarks have the following decomposition 

h = aiP\ + PiP2 + k it x, i = 1,2. (4) 

Taking into account the on-shcllness of the produced quarks, the above Sudakov parameters can be 
expressed in terms of rapidities, transverse momenta and heavy quark masses M, i.e. 



M 2 + k 2 I M 2 + k? 
-e v \ 



A = Y ^ e "" i ' i = 1 » 2 ' (5) 

Here 771 (772) is the rapidity of the produced heavy quark (anti-quark) and kf = —kf ± are the corre- 
sponding Euclidean squared transverse momenta. 
Making use of the definitions 

si = (pi + q 2 ) 2 = x 2 a, s 2 = (p 2 + qi) 2 = x\ s, (6) 

which correspond to the center of mass energies of the upper and lower subamplitudes in Fig. [TJ 
respectively, we can write the following expression for the differential cross-section of heavy quark 
production: 



d 6 a f 1 , f 1 , f d 2 qi f d 2 q 2 



= / dx\ I dx 



drndmd 2 kid 2 k 2 J 'Jo *J (2tt) 3 J (2tt) 3 

l r RR^Qq(qi:q2;ki,k 2 ,z)| 2 



d 2 q a ^^(qa) , / si 



■> f[- — ,q a ,qi 



X 



2 2 

qfq2 



2tt <\l V s o,i 
d 2 q b <J> B (q 6 ) / s 2 



/ — ,q2,qb 



2tt q 2 * V s o,2 : 
x (2^) 4 <5( 2 )( qi + q 2 - ki - k 2 ) S( Xl - ai - a 2 ) 6(z 2 - ft - A). (7) 

In this expression $^4 and <I>b denote the hadron impact factors, which are responsible for the coupling 
of the reggeised gluon to the proton A and B, respectively. r RR ^QQ indicates the high energy effective 
vertex coupling the two reggeised gluons to the heavy quark-antiquark pair with 

xi ^/kf + M 2 + Vk| + M 2 e^-vi 

being the fraction of the longitudinal momentum of the upper reggeised gluon along pi, carried by 
the heavy quark. / denotes the BFKL gluon Green function with the following Mellin transform / w : 

/(— ,q ,qi)= / ( — ) / w (q a! qi), (9) 



so,i J Jc 2tt« \ s o, 

where the contour of integration C lies parallel to the imaginary axis and to the right of all the 
singularities in f u . The resummation of high energy logarithms is achieved by iterating the BFKL 
integral equation for f u : 

w/ w (q OI qi) = <^ (2) (qa - qi) + J rf 2 q^BFKL(qa,q)/^(q,qi)- (10) 

In a general case where both the produced heavy quarks and the impact factors would provide a 
similar hard scale (i.e. if the protons were replaced by highly virtual photons, or jets with a high p t 
were tagged in the forward/backward regions), a good choice for the energy scales so,j would be given 
by so,i = |qo|VS and s , 2 = VS|q b |, where 

S = x x x 2 s = s+(ki+k 2 ) 2 , s = {k x +k 2 ) 2 . (11) 
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s reads for the squared center of mass energy of the partonic process g*g* — > QQ. Such a choice 
naturally introduces the rapidities rj^ and r/g of the emitted particles with momenta and since 

— = e fe-'7QQ) w ) ( -^-\ = e&QQ-vs)" ; (12) 

s 0,l/ V s 0,2/ 

with the rapidity of the heavy quark system being given by 

1, ai+a 2 , . 

^ = 2 k ftTft' (13) 

If the BFKL Green function, the impact factors and the production vertex were known exactly at NLO, 
then the precise choice of the energy scales So,i would turn out to be irrelevant since any dependence 
of the cross section on this scale would cancel at the same NLO accuracy. However, even in that case, 
when the kernel is exponentiated there is a residual dependence on so,i which would correspond to 
NNLO and higher terms. A natural choice for so,i is then that which reduces the size of those higher 
orders corrections to the minimum for a given observable. 

In the case of interest for us in this letter there exists a hierarchy of scales with a large difference 
between the only hard scale provided by the invariant mass of the heavy quark pair system and 
the large transverse size of the incoming hadrons. Here the previous symmetric choice of scales is 
not appropriate as the scale of the heavy quark anti-quark system E is significantly larger than the 
transverse scales and associated to the scattered protons. A more natural choice for so,, is given 
by E alone, i.e. 

'— V - *T", (— 

This choice of the energy scale is common in deep inelastic scattering and leads to the concept of the 
unintegrated gluon density in a hadron. This represents the probability of resolving an off-shell gluon 
carrying a longitudinal momentum fraction x off the incoming hadron, together with a transverse 
momentum kx- 

As it is well-known, any choice of energy scale only matters at next-to-leading and higher orders 
since the LO approach is scale invariant. The LO unintegrated gluon density g LO is defined as 



(14) 



,LO 



= r !s*££0/lo ( k)l with /L o (X;qjk) = i ^^ /L o (qk)) (15) 

J 2tt q 2 J c 2m 



where corresponds to the solution of the LO BFKL equation with kernel -RTbfkl = -Kbfkl- 
Contrary to the LO case, the next-to-leading order BFKL evolution is sensitive to changes in the 
energy scales so,i- As it was pointed out in [21], any shift of scales can be absorbed in the kernel, 
impact factors, and central production vertex. With the choice of energy scale as in Eq. (TTJJ the NLO 
impact factors are modified by an extra logarithmic term of the form 

^ LO (q) = *NLO(q) - i [ J^O^q) hX (16) 



J2 br lVL \ 5 17 — q2 • 

The NLO kernel receives two additional contributions, corresponding to the incoming and outgoing 
reggeised gluons: 



^BFK L (l ai h) =K^° h (l a , h)~\J d 2 \K^{\ a , 1) i^° KL (l, h) In jj. (17) 



The NLO QQ production vertex also gets two types of corrections, corresponding to the two different 
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evolution chains originating from the hadrons A and B: 

|rRR°QQ(qi,q 2 ;ki,k 2 ,z)| 2 = |r RR ° QQ ( qi ,q 2 ; kl ,k 2 ,z)| 2 

- f / ^l^ Q Q(qi,l;k l ,k2,*)| 2 «-^KL(I,c tl )ln^ 



I) 2 ' 

With these modifications, the NLO unintegrated gluon density is defined as follows 



(18) 



g^°(x,k) = [ ±^^mf {x , q ,k), />,q,k) = /^-^ (q , k ), (19) 

J 2tt J c 2m 

where / w obeys the modified NLO BFKL equation 

w/a,(q*,qi) = <5 (2) (qa - qi) + / d 2 q^BFKL(q Q , q)/ w (q, qi) (20) 
with a NLO kernel which we will discuss in Section 3: 

/W L (q Q , q) = K™ Kh (q a , q) + ^g L (q a , q). (21) 

Using these definitions, the differential cross section in Eq. {7J at NLO accuracy is given by the 
expression 

d6<7 f 1 j f 1 j /" rf2 qi f rf2 q2 NLO/ \ NLO/ ^ 

dxi / cir 2 / j^—^t / t^-tt 5 (^i,qi)3 (^2^2) 



d m d m d?\i l( P\i 2 J L J Z J (2tt) 3 7 (2tt) 



r RR^QQ(qi.q2;ki,k 2 ,z)| 



2 



x 1 rr^qq 2 2' <5 (2) (qi + q 2 - ki - k 2 ) S( X1 - ai -a 2 ) 6(x 2 -fr- fa) (22) 

qiq 2 

rf2 qi nlo/ x nlo/ , a , J r RR^QQ.(qi> k i + k 2 - q x ; ^ , k 2 , z) | 2 



/ ^- 9NL ° (ai + "2, qi)3 NLO (/3i + 02, ki + k 2 - qi 



(2tt) 6 qf(ki +k 2 - qi) 2 

This expression can be interpreted as the convolution of the unintegrated gluon densities with a 
partonic differential cross section i.e. 

d 6 a f 1 dxi f 1 dx 2 f d 2 qi f d 2 q 2 NLO , , nlo/ _ \ 



>(x 1 ,q 1 )g N ^(x 2 , q2 )- 



drjidrj 2 d 2 'kid 2 'k 2 J X\ J x 2 J 2ir 2 J 2ir 2 ' ' dr)idri 2 d 2 \iid 2 k 2 

(23) 



where 



d 6 * - il^^ 2 ^- fc2 )| 2 ^ 1 ^^( 2 -) 45(4) fe +32 - fer - fa). (24) 
Here 2E is the flux factor and 

S 2 

l-4(9i, 92; fci,fe)| 2 = ^- T |r RR ^ Q q(qi,q 2 ;ki,k 2 ,z)| 2 (25) 

qiq 2 

is the squared matrix element for the production of a heavy QQ pair from the fusion of two transversely 
polarised reggeised gluons. This means that their polarisations are chosen to satisfy 

E e (A)fe)<A)fe) = ^r-, with i= 1,2, (26) 

\ q« 



5 



and can be related up to the overall factor £ 2 /q 2 q 2 in Eq. (|25|) to the usual longitudinally polarised 
reggeised gluons by means of a Ward-identity for the i-channel gluons. 

At present the NLO corrections to the heavy quark production vertex r RR _,QQ are not available 
and only the LO vertex is known . This LO result can be written in the following form 



l r RR^QQ(qi'q2;ki,k 2 ,z)| 2 = 5 4 f ^A 1 (q 1 ,q 2 ;k 1 ,k 2 ,z) + ^A 2 (qi,q 2 ;ki,k 2 ,z) ) , (27) 



with 



Ai(qi,q 2 ;ki,k 2 ,z) =-— 2' 



sE \(^A/2 u-M 2 )\ z l-z 

[(qi - kQ 2 + M 2 ] [(q x - k 2 ) 2 + M 2 ] - (k 2 + M 2 ) (k 2 + M 2 ) \ 2 

(t-M 2 )(u-M 2 ) 
(q 1 -k 2 ) 2 +Af 2 - 7 ^(k 2 +M 2 ) E(M 2 ) 
u-M 2 + s 

< / (q 1 -k 1 ) 2 +M 2 -^(k 2 + Af 2 ) E(M 2 ) 
X V i-M 2 s 

and 

[(qi - kx) 2 + M 2 ][( qi - k 2 ) 2 + M 2 } - (k 2 + M 2 )(k 2 + M 2 ^ ^ 2 



(28) 



A2(qi,q2;ki,k 2 ,z) = 



(i - M 2 )(u - M 2 ) 



q?qi 



(t — M 2 ){u — M 2 ) 



(29) 



To write down Eqs. (|28I29[) we defined, apart from the variables introduced in Eqs. (|8I11[) and the 
partonic Mandelstam invariants 

i = ( qi - h) 2 = -^04 + m 2 ) - ( qi - ko 2 , (30) 

fi = («i - k 2 ) 2 = -^(k 2 + M 2 ) - ( qi - k 2 ) 2 , (31) 

1 — z 

the following set of transverse momenta 

A = ki+k 2 , A = (1 - *)ki - *k 2 , (32) 

which allow us to express Eq. (JTTJ) as 

A 2 4- M 2 

E = s + A 2 = \ + A 2 . (33) 

z(l - z) 

Finally, we have also used 

E{M 2 ) = 2(2z - l)q 2 + 2qi • A + \~ 2Z (A 2 + Af 2 ) - [(2z - 1) A 2 + 2A • A] (34) 

The explicit form of the vertex in Eq. (|27p , keeping all the information on the outgoing QQ system, 
will permit a comprehensive study of differential distributions in exclusive observables. For this we 
will also need to keep track of the multiple soft emission stemming from the gluon evolution. How to 
achieve this task is discussed in the following Section. 
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3 Multiparticle production and the unintegrated gluon den- 
sity 

The NLO unintegrated gluon densities in Eq. (fl9|) , which enter the differential cross section of Eq. (f22|) , 
require both the NLO BFKL gluon Green function and the proton impact factor <I>p(q). The latter is 
of non-perturbative origin, it can only be modelled and has to be extracted from the data. A possible 
simple choice for a model of the proton impact factor would be 

^-(^) A - < 35 > 

Here A is a positive free parameter, while A is a momentum scale of the order of Aqcd- A more 
sophisticated alternative to Eq. (|35[) has been presented in Ref. [22] where it was proposed to expand 
the proton impact factor over a set of orthogonal conformal invariant eigenfunctions. 

The second building block of the unintegrated gluon densities is given by the NLO BFKL gluon 
Green function. An alternative formulation to the usual treatment in Mellin space was proposed in 
Ref. [23j . This form of solving the equation by iteration in momentum space has the advantage of 
dealing exactly with running coupling effects and incorporates the full azimuthal angle dependence 
of the soft multiparticle emission in multi-regge kinematics associated to the BFKL evolution. For a 
complete analysis of exclusive properties of multigluon final states associated to the production of a 
heavy QQ pair we consider this to be the most convenient of the available methods of analysis of the 
BFKL Green function at NLO. By means of a phase space slicing parameter A the virtual and real 
contributions are treated separately. 

As it has been previously explained, the BFKL kernel receives in the case of NLO unintegrated 
gluon densities an additional contribution, Eq. (|17[) . due to the choice of the energy scales So,i. This 
affects the real emission contribution to the kernel but not the gluon Regge trajectory, u\, which in 
this physical regularisation can be written as 

WA(q) = -f(|q|A)ln^+^|c(3) with = a s + f (| - ^ + |^ - <£- lu^j , (36) 

where a s = a s (/i)A c /7r and (3q = (HJV C — 2n/)/3. [i is the renormalisation scale in the MS scheme. A 
can be understood as an effective gluon mass or as a lower cut-off for the transverse momenta of the 
emitted gluons. 

The NLO real emission kernel K\ , which accounts for the emission of gluons or massless quarks 
in quasi- multi-regge kinematics, is given by the following sum 

Kl ea \l a , l a + 1) = -U( l2 M 12 - A 2 ) + K™\l a , la + 1) + K coll (l a , L + 1), (37) 
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where [20] 



111 ' n*J { iei 2 i 2 



l 2 ' 1? ' W 2 \l 1 



(1 2 +1 2 ) 2 (21 2 1 2 - Sl^-Sl 4 ) 



16141^ 



(7.T 



1 



l 2 +x 2 l 2 



In 



1 + x 



2(1 2 "1 2 ) 



(l a -l b ) 2 (l a + l b ) 2 



1-Z 



(l^-l 2 ) 2 



(l a -l 6 ) 2 (l a + l 6 )< 




) t 



! ;^(i b -zi Q ) 2ii± i 2 



In 



(38) 



and 



Kcou(W&) = " ^ y d 2 l^K L (la,l)^BFK L (l,^) ^ p. 



(39) 



Following Ref. J53j, this representation of the NLO BFKL kernel can be now used to solve iteratively 
the integral equation. The explicit solution for the gluon Green function then reads 

{oo n r i- i— 1 i 

*< 2 >( q - k) + J2 n J [^ cal (q + E 1 ^ + E w 



n=l i—1 1 

dXj -u x (q+J2' j = 1 ljO+^Afq+E-;! lj) 

3*2 



-5 (2) (q + E^- k )h ( 4 °) 



with xq = x. Note that this representation can be now implemented in a Monte Carlo event generator 
where all the information about each of the emitted particles is recorded. At NLO each iteration 
of the kernel, or each of the terms in the sum of Eq. (|4"0")) . corresponds to one or two emissions well 
separated in rapidity from previous and subsequent clusters of particles. Inserting this function in the 
formula for the differential distributions will generate our exclusive observables. 

In the real emission kernel of Eq. (|3T|) . the two terms explicitly written in Eqs. (|38l39p do not carry, 
apart from an overall <5 2 (/z 2 ) factor, any renormalisation scale dependence. This is different from the 
remaining part of the real emission kernel and the gluon trajectory which contain the function £ of 
Eq. (|36p. This can be written as 



Z(X) = Mm 2 ) i 



a,(/i 2 ) a . x 

— — — /3 m — 



a s (fJ, )S I , with 5 = 



1 

12 



.(3o 



(41) 



In this expression, the logarithmic term can be absorbed into a redefinition of the running of the 
coupling which corresponds to the replacement of a s {p?) by a s (X). The remaining, non-logarithmic 
term, can be identified as a common factor which appears when dealing with resummations of soft 
gluons [24;. The term a s (1 + a s S) is proportional to the two-loop cusp anomalous dimension. Gener- 
ally, the appearance of this term offers the possibility to change from the MS renormalisation scheme 
to the Gluon- Bremsstrahlung (GB) scheme. Such a change corresponds to a shift of the Landau pole 
Agb = Ajjg-exp S?j^-. Stability under this change of scheme offers a good tool to test our theoretical 
predictions. 
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4 Conclusions 



In this letter we set the theoretical framework for a study of heavy flavor production in central 
regions of rapidity in hadron-hadron collisions using kx factorisation at NLO. While the heavy flavor 
production vertex is kept at LO, the unintegrated gluon density is treated by taking into account 
the full NLO corrections. The latter contain both the full NLO BFKL evolution and the further 
corrections which arise due to the asymmetric choice of energy scales, inherent to hadronic cross- 
sections, and which can be understood as the onset of collincar evolution from the soft hadrons to the 
hard production vertex. 

The NLO BFKL Green function which, convoluted with the (non-perturbative) proton impact 
factor, forms the NLO unintegrated gluon density has been presented in an iterative way which allows 
for a numerical evaluation using Monte-Carlo integration techniques. In future publications we will 
present results on this numerical implementation, together with fits of our unintegrated gluon density 
to deep inelastic data from HERA and predictions for heavy quark pair production at the Large 
Hadron Collider at CERN. 
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